A NEW ALGORITHM FOR THE RECURSION OF HYPERGEOMETRIC 
MULTISUMS WITH IMPROVED UNIVERSAL DENOMINATOR 



STAVROS GAROUFALIDIS AND XINYU SUN 

Abstract. The purpose of the paper is to introduce two new algorithms. The first one computes a linear 
recursion for proper hypergeometric multisums, by treating one summation variable at a time, and provides 
rational certificates along the way. A key part in the search of a linear recursion is an improved universal 
denominator algorithm that constructs all rational solutions x{n) of the equation 

-——x{n + m) H h ——x{n) = c(n), 

bm{n) bo{n) 

where ai{n),bi{n), c{n) arc polynomials. Our algorithm improves Abramov's universal denominator. 
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1. Introduction 

1.1. History. The paper introduces a new algorithm to find linear recursions (with coefficients polynomials 
in n) for multidimensional sums of the form 
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(1) 

keD 

where _D C Z*" and the summand /(n, k) is a proper hypergeometric term in the variables (n, k). By proper 
hypergeometric term (abbreviated by term) f{m) in the variables m — (mi, . . . , mg) we mean an expression 
of the form 



J 

(2) f{m)^P{m)l[A,im)\^^ 

where P{rn) is a polynomial in ?ti and Aj{m) = cijirrij is a linear form in m with integer coefficients 

and = ±f for f < J < J. Throughout this paper, f{n, k) will denote a proper hypergeometric term. 

As observed by Zeilbcrger |Ze| . and further explained in |WZ| . Sister Celine's method |Fas| can be used 
to prove the existence of linear recursions of S{n) in a constructive way. A faster algorithm was constructed 
by Zeilberger (also known as creative telescoping [PWZj ) . which employed Gosper's indefinite summation 
algorithm |Gos| . Creative telescoping is faster than Sister Celine's method, and often returns the optimal 
(i.e., minimal order) recursions. However, due to the nature of Gosper's algorithm, Zeilberger's method only 
works for single sums, i.e., when r = f in ([T]). 

Wegschaider in |Weg| improved Sister Celine's algorithm for multisums; Zeilberger has a program EKHAD 
for creative telescoping, while Paule and Schorn ^PSh^ implemented it in Mathematica; Schneider created 
a package called Sigma, the framework of which was explained in [Sch| : Apagodu and Zeilberger |AZ] 
generalized creative telescoping to multi-variable context which resulted in another fast algorithm. 

1.2. What is multivariable creative telescoping? Multivariable creative telescoping for S{n) is the 
problem of finding a natural number J S N, and rational functions aj{n) £ Q(n). for I < j < J and rational 
functions Ci{n, k) S Q(n, /c) for 1 < i < r so that 

J r 

(3) ^a, (n)7VV(n,fc) = ^(A^, - f )(Q(n, A:)/(n, fc)), 

j=0 i=l 

where N ., Ki, n, ki are operators that act on functions ,f{n, k) as follows: 

[Nf)[n,k) = /(n+f,fc), {nf){n,k) = nf{n,k), 

(4) {KJ){n,ki,...,kr) = f{n,ki,...,ki_i,k, + l,ki+i,...,kr), 

{kj){n,ki, . . . ,kr) = hf[n,ki,...,kr). 

Note that the operators iV, n, Ki, ki commute except in the following instance 

(5) Nn = n + 1, Kiki = k^ + 1. 

In Equation ([3]), the rational functions Ci{n,k) for 1 < i < r are called the certificates and the operator 
X]j=o o,j{'n')N^ is called the recursion for the sum S{n). Given Equation ([3]), we can sum over k to obtain an 
inhomogeneous linear recursion for S'(n), whose inhomogeneous part consists of the contribution from the 
boundary terms. 

All known algorithms of creative telescoping convert ([3]) to a system of linear equations with coefficients in 
the field Q(n, k). This is possible since dividing both sides of ^ by /(n, k) and using the fact that /(n, k) is 
proper hypergeometric, it follows that the ratios Nf{ji, k)/ f(n, k) G Q(?t., k) and Kif{n, k)/ f{n, k) G Q{n, k) 
are rational functions. The number of unknowns and equations directly affect the performance of the above 
mentioned algorithms. 
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1.3. Abramov's universal denominator algorithm. A key part of our search for a linear recursion of 
hypergeometric multisums is an improved universal denominator algorithm that finds all rational solutions 
x{n) £ Q{n) to a linear difference equation 



amin) . , ^ , , ao{n) 

- — r^xin + m) + ■ ■ ■ + xin) = c(nj, 

bm[n) bo(n) 



where ai(n), bi(n), c{n) are polynomials. The idea is to correctly predict the denominator u(x) of xin) (also 
known as the universal denominator) , so that the problem can be reduced into finding a polynomial solution 
to a linear difference equation. In |Abj Abramov developed a universal denominator algorithm. In this 
paper, we develop a new algorithm that improves Abramov's algorithm by possibly reducing the number of 
factors in the universal denominator. The new algorithm is used repeatedly to convert the problem of finding 
recursions of multivariate hypergeometric sums into the problem of solving system of linear equations. And 
fewer factors in the universal denominator implies fewer numbers of variables and fewer equations in the 
system. 

1.4. Acknowledgement. The authors wish to thank D. Zeilberger for a careful reading of an earlier version 
of the paper and for detailed suggestions and comments. 

2. Two ALGORITHMS 

2.1. A new algorithm for the recursion of hypergeometric multisums. To describe our algorithm 
for the recursion of multisums, let us introduce some useful notation. 

Definition 2.1. Fix a term f{n, k) where fc = (fci, . . . , kr) and 1 < i, j < r. We say that two operators P 
and Q in the variables n,ki, N and Ki are /-equivalent modulo Ki ^ 1, . . . , Kj — 1, and write 

(6) P^sQ mod [K, - 1, K,+i - 1, . . . , - 1), 
if there exist rational functions k) for i < s < j so that 

(7) (P - Q)/(n, k) = ^(A', - l){bs{n, k)f{n, k)). 

s—i 

If i > j, the right-hand side of the last equation is 0. 

Our algorithm will construct operators RECi for < i < r and RECj^ for < i < j < r of the following 
form 

di do 

RECi = ^a,,,(n,fci,...,fc,)^,', « ^ 0; RECo = J]] "o.^N A^; 

s=0 s=0 

dj,i djo 

RECj,i = A, + J2 aj.iAn, ki,..., kj)K], i ^ 0; RECj,o = N + J2 ajfiA^, ki,..., kj)K], 

s=0 s=0 

that satisfy 

(8) RECi mod {Ki+i - 1, . . . , A^ - 1) RECj,i =/ mod (A^+i - 1, . . . , A^ - 1). 

We will call such operators RECi,RECj.i f -compatible. 

Observe that RECo is exactly Equation jS]). Here are the steps for the algorithm. 
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INPUT: A proper hypcrgcometric term /(n, ki, . . . , kr). 

OUTPUT: A recursion X]i=o ^iC'^)-^' ^^'^ ^ set of certificates Ci(n, fci, . . . , kr) that satisfy ([3]). 



Step 1. 
Step 2. 

Step 3. 
Step 4. 
Step 5. 
Step 6. 
Step 7. 



1 < j < r - I 



Set I := r, ko = n and /io 
Set RECr^o :^N-^, RECr,i := / 
and RECr := A"^ - 

Construct RECr-i using Proposition 14.41 

If / = 1, print RECo and stop; otherwise, continue; 

Construct RECi-i^ for < i < Z — 2 using Proposition 

Construct RECi_2 using Proposition 14.41 

Set 1 = 1 — 1, and go to Step 4. 



There is some similarity between our algorithm and results of Schneider [Sch| : we do believe however the 
underlying algorithm to obtain the certificates is different from Schneider's program Sigma, although he did 
employ some version of Abramov's algorithm. 

The subtle part of the above algorithm are steps 5 and 6 which compute the proper denominators for the 
certificates that appear in Equations ^ . This is done using Propositions 14.31 and 14.41 which follow from 
Theorem 15. 1[ which are implemented in our improved denominator algorithm of Section [ 

Example 2.2. When r 
following order: 



3 the algorithm computes RECj for < j < 3 and RECj j for < j < i < 3 in the 



REC3 0, REC3 1, REC3 2, REC3 — > REC2 > REC2 Oj REC2 1 — > RECi — > RECi — ^ RECq 



A Maple implementation of the above algorithm is available at [GS2j . A Mathematica implementation 
will be developed later. A g-version of the above algorithm is possible and will also be developed later. 

2.2. An improved universal denominator algorithm. In this section we present our universal denom- 
inator algorithm. Let K denote a fixed field, which in applications it is the field of rational functions with 
rational coefhcients in a finite set of variables. 



INPUT: An equation with rational coefficients x{s + m) + • • • + x{s) = c(s), 

where ai[s),bi{s), c{s) 6 K[s] are polynomials. 
OUTPUT: A rational solution x{s) = ^^^^ to the equation, where P(s), Q(s), R{s) £ K[s] 



Step 1. 
Step 2. 
Step 3. 
Step 4. 
Step 5. 
Step 6. 
Step 7. 
Step 8. 

Step 9. 

Step 10. 
Step 11. 
Step 12. 



Set P{s) = 1; 

Set (t(s) = 1cm (6o(s), , 

Set Q{s) = gcd(ro(s), . 



^{s)), and Ti(s) = a,(s - j)^^^, 



< i < m; 



Find the largest possible nonnegative integer £ such that gcd(^^, ^"^^['J)^ 
If such an £ does not exists, continue to Step 9; 
Otherwise, set Q{s) = Q{s) YII^q (I){s + i); 
Set To(s) = and r„,.(s) 
Go to Step 4; 

Set R{s) = lcmo<,<,,„ | g,d(6,(«-,)Q(s),a.(s-z)di'J)?(itn,^, fc.(''-«)Q(^+J-0) 

If = 1, STOP; 

Otherwise, set P{s) = P{s)R{s); 

Go to Step 9. 



(s) ^ 1; 



2.3. Plan of the proof. The structure of the paper is as follows. In Section [31 we explain the usage of the 
Maple program. We present a few examples and compare the results and performance against the programs 
discussed above. In Section O we introduce the terminology used in the paper, and present the general 
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structure of the method as a sequence of steps. In Section HI we prove the vahdity of each step of the 
structure, and also explain the method in detail. In Section [51 we prove a new algorithm that generates 
universal denominators with possibly less factors than those generated by Abramov's algorithm, that also 
partially predict the numerators for rational solutions to linear difference equations. 

3. Use of the program and examples 

Example 3.1. Define 

g{n, k) = ( " 

We will prove that f fPWZl Page 33] and [Sch]) 

^ f{n,ki,k2) = ^g{n,k). 

After running the program, both sides of the above equation are annihilated by the operator 

(n + 2fN^ - 2{2n + 3){3n^ + 9n + 7)N - 4(4:71 + 5)(4n + 3)(n + 1). 

Since they have the same initial conditions for n = 0, 1, the two sides agree for all natural numbers n. 
Please see |GS2| for the syntax of input and output. 

ExaTnple 3.2. Define 

n \ ( n\ fn + ki \ f n + k2\ /2r7, — fci — k2 



2 



yklj \fc2/ \ / V ^2 / V 

Please see }GS2| for complete information. The recursion for the multisum '^^ki k2 '^i' ^2) is of degree 4. 
ExaTnple 3.3. Define 

/,„, M . [I) ) ) (» ;;■■) (" (» (- - - - 

Please see |GS2] for complete information. The recursion for the multisum ^^.^ ^.^ ^.^ /(n, fci, A:2, /ca) is of 
degree 4. 

4. Proof of the multisum algorithm 

4.1. Two Lemmas. We fix a term f{n, k) where k = (fci, . . . , kr), and consider a fixed variable ky and the 
corresponding operator Ky. The moduli are always . . . , ky), which we suppress for simplicity. 

Lemma 4.1. //X]i=o ^i: ■ • ■ j ky)Kl =/ aTid N + X]i=o ^^ii^T^i: ■ ■ ■ ^ ky)Kl =/ 0, then for any integer 
m and rational functions {ai{n, fci, . . . , A:„)}o<i<m; there exist rational functions {/3j(n, fci, . . . , fct,)}o<j<7-i 
so that 

m 7-1 

(9) ^ a,(n, fci, . . . , ky)N' ^ /?, (n, fci, . . . , fc-„)i^^. 

i=0 j=Q 

Furthermore, Equation ([9]) is a linear system of equations with unknowns {/3j(ri, fci, . . . , ky)} and coefficients 
in the field Q(n, fci, . . . , fct,). 

Proof. Since the operators are linear over the field Q{n, fci, ... , fc„), we only need to show the result for N"^ 
for m > 1 by induction. 
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The conclusion is true for m — 1. Suppose it is true for m — 1, i.e., N"' ^ =/ Sj=o 7i('^i ^i' ■ ■ ■ ' ^v)Ki 
for some {7^}. Then we have 

/-I i-i 
iV™ N{N"^-^)^fNC£l,{n,ku...,k,)Ki)^fY,j,{n + l,k,,...,K)KlN 

j=o j=0 
i-i / i-i \ i-i 

=/ ^lj{n + l,ki,. . . ,k^)Kl I ~^ai{n,ki,...,ky)Kl j =f ^ fci, . . . , fc„)7^*, 

j=0 \ i=Q ) i=0 

for some rational Pi{n, ki, . . . ,ky). The last equation is because the order of recursion in k^ satisfied by / is 
at most /. 

Since the reduction of N"^ does not depend on {ai{n, fci, . . . , ky)}, it follows that {/3j(n, ki, . . . , ky)} are 
linear functions of {ai{n, fci, . . . , ky)}. □ 



Lemma 4.2. Given + YTi^a ai{n., ki, . . . , ky)Kl =/ and 

p—X p— 1 

-^ap_i_j(n, fci,...,fc„+j)6p_i(n, fci,...,fc„ + l+j)-6p_i(n, fci,...,A:„) = ^ Cp_i_j(n, fci, . . . , fc„ + j), 

where {ci(n, fci, . . . , A:«)}o<i<p-i a"''^ bp-i(n, ki, . . . , ky) are rational functions. Define, for Q < i < p — 1, 

bi{n,ki, . . . ,ky) = bp^i{n,ki, . . . ,ky - p+1 + i) + ^ aj+j(ri, fci, . . . , fc„ - j)6p_i(n, fci, . . . , fc^, - j + 1) 

p — 2— 1 

+ X! c.i+j{n,ki,...,ky - j). 
i=i 

T/ien 

p— 1 p— 1 

c,(ri,, fci, ... , fc„)fci =/ (ky - 1) ^ bi{n, fci, ... , fci,)fc^ 

1=0 i=0 

Proof. From the definition of bi{n, fci, . . . , fc„), it is easy to check that 



6i_i(n, fci, . . . , fc„ + 1) = bi{n, fci, ... , fct,) + ai{n, fci, ... , ky)bp^i{n, fci, ... , fci, + 1) + Ci(n, fci, ... , fc,;), 
boin,ki, ...,ky) = -6p_i(n, fci, . . . , fc-„ + l)ao(n, fci, . . . , fc^,) - Co(n, fci, . . . , fc,;). 



It follows that 
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p-1 p-2 

{K, - 1) 51 ^«("' '^i' ■ • ■ ' =f ^P-i(«> + l)KP + J2 Hn. h,...,k, + l)Kl+^ 

p-1 

-^6i(n, ki,...,ky)Kl 

i=0 

p-1 

=f -^6p_i(n, fci, . . . , fey + l)ai{n, fci, . . . ,ky)Kl 

i=0 

p-1 

+ ^ (6i_i(rt, ki,...,ky + l)- bi{n, fci, . . .,ky))Kl - bo{n, fci, . . . ,fc„) 

i=l 

p-1 

=f -^6p_i(n, fci, . . . ,fc„ + l)aj(n, fci, . . . ,ky)Kl 

i=0 
p-1 

+ '^{ai{n,ki, ky)bp-i{n,ki, ...,ky + l) + Ci{n, fci, . . . ,fc„)) 
1=1 

+6p_i(n, fci, . . . , fci, + l)ao(n, fci, ... , fc„) + Co(n, fci, . . . , fci,) 

p-1 

=/ ^Ci(n, fci,...,fc„)X^. 

i=0 

□ 

Lemma 221 also appeared in jSchj in a different form. It is included here for completeness of the proofs. 



4.2. Two propositions for the algorithm. In this section we state and prove Propositions 14.31 and 
which are used in our algorithm. Fix a term f{n,k) where fc = (fci,...,fcr). Recall we set ko ~ n and 
Ko = N. 

Proposition 4.3. Let 1 < v < r. Given f -compatible operators RECv+i, RECv, RECv+i.u, R'ECv+i,v 

for < u < V < r, it is possible to construct an f -compatible operator RECv,u for < u < v in Step 5. 

Proposition 4.4. Let \ < v < r . Given f -compatible operators RECv, RECv, v-i for < u < v , it is 

possible to construct f -compatible operator RECv-i in Steps 3 and 6. 

4.3. Proof of Proposition [4T3l Let 

(10) 

,7-1 

RECv+i := A";,+i + ^ ai(fco, fci, ... , fci,+i)i^^+i and RECv+i =/ mod {L<y+2 - I, . . . , Kr ~ I). 

j=0 

We can always divide the operator by the leading coefficient if it is not 1, since it does not involve variables 
fc^-i_2, . . . ,kr. Let us look for 

7-1 

(l]9,ECv,u Ku + Mko, ki,..., k^)Ki and RECv,u mod (/v„+i - 1, . . . , - 1), 

1=0 

for some rational functions ^^(fco, fci, . . . , fc«). To prove the existence of RECv.u, borrow the idea in the 
proof of [WZ| by solving 



(12) 



t7„ U„ 



Aru + X! "' X! '^i.--t^{'^o,ki, . . . ,ky)Y\_Ki 

i„=0 im.=0 \ l=v 



/(fco,fci, . . . ,fcr) = 0, 
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with ai^...i^{kQ, ki, . . . , k^) being the unknown rational functions. Divide both sides by the hypergeometric 
function /(fco, . . . , kr) to obtain an equation of rational functions. By comparing the coefficients of the powers 
of fci,+i, . . . , kr^ we can set up a system of linear equations over the field Q(fco, fci, . . . , ky), whose unknowns 
are ai^...i^^ (feg, fci, . . . , ky). The number of unknowns is ni"i)(^i + 1); while the number of equations, which 
equals the degree of the numerator in Equation IT^ . is proportional to (n"=t; ^i) {j2tLy jj;") • It follows that 
when Uy, . . . , Um are large enough, we have more unknowns than equations in the system, which guarantees 
a nontrivial solution. Replacing ifu+i, . . . ^Kr in Equation (jl2p with 1, we get a solution to Equation (jlip . 
The maximum power of J — 1 on Ky is ensured by the existence of a recursion of order /. The readers may 
also compare with |PWZi Theorem 4.4.1] or [AZ[ Theorem MZ] for a detailed discussion on the method in 
similar cases. With the proof of existence completed, we can introduce a new method to find the functions 
{4>i} and {bj}. 

Reduce + YhZI 4>i{ko,ki, . . . , ky)Kl into Yl,iZo Ci{ko, ki, . . . , ky+i)Kl^-^ for some rational 
Ci(fco, fci, . . . , ky^i), using Lemma |4. II below. This implies that 

^ c^{ko, ki, . . . , ky+i)Ky_f_i =f {Ky+i - 1) I ^ b,{ko, ki, . . . , ky+i)Kl^i \ . 

i=0 \i=0 I 

Since the coefficient of Ky^^ is 1 in pO|) . it follows from Lemma 14.21 below that we only need to find 
fe,7-i(fcoj fcii • ■ • J ^u+i) such that 

,7-1 

- X] "./-i-ie^o, ^1, ■ • ■ , K+\ + j)bj-i{ko, fci, ... , fct,+i + 1 + j) 



,7-1 



(13) = Xcj_i_j(fco,fci,...,fct,+i 

In the equation. {ai(fco, fci, ... , fc,u+i)}o<i<,7-i a-re known; bj_i is a rational function over fco, . . . , fc^+i over 
the field Q(fco, fci, ■ • ■ , ky); and {ci}o<i<,/-i are linear over {0j}o<j<7-i- So the right-hand side can be written 

as ^j=o ^''^''°y^l'^'^',^^^^^^^^ , with polynomials 0j(fco, fci, . . . , fc„) unknown; and C/j(fco, fci, . . . , fc-i,+i) 

and V{kQ, fci, ... , fc^+i) known. Multiply both sides of Equation (|13p by ^(fco, fci, ... , fci;+i) to obtain 

j-i 



-aj_i_j(fco,fci, . . . ,fc„+i +j)bj-i(ko, fci, . . . ,fc„+i + 1 + j)T^(fco,fci, . . • ,fc«+i) 

i-i 

(14) = X C/j(fc-o,fci, . . .,ky+i)(j)j{ko,ki,. . . ,ky). 



3=0 

In the above equation, consider 6j-i(fc'o, . . . , fc^j+i) G Q(fco, ■ • ■ , kv){ky+i), and apply Theorem [ST] to the field 
K = Q(fco, ■ • ■ , fcu) and the variable s — ky+i. It follows that we can write 

i?(fco,fci,...,fct,+i)P(fco,fcl,-.-,fci,+l) 

6j_l(fco,fcl,...,fc„+l) = —T- ; ; r , 

y (fco, fci, . . . , fc'u+l) 

with polynomials i?(fco, fci, ... , fct,+i) G Q[ko, . . . , fct,+i] and Q{kQ, fci, ... , fct,+i) £ Q[fc"o, • • ■ , fct>+i] known, and 
P(fco, fci, ... , fc^+i) G Q[fco, ■ • ■ , fcu+i] unknown. By multiplying both sides by the common denominator of the 
left-hand side, and comparing the degree of fc„+i, we can determine the degree of fcu+i in P(fco, fci, . . . , fctj+i), 
say, L. By writing P(fco, fci, ... , fct,+i) as J2i=o V'i(fcoj fci, ■ • ■ , fcii)fci+i, plugging it back into Equation (|14p . 
and comparing the coefficients of powers of fc^+i, we can set up a system of linear equations with {t/)j}o<j</-i 
and {ipi}o<i<L as unknowns. The system is guaranteed to have a nontrivial solution because of the existence 
of the recursion. □ 
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4.4. Proof of Proposition 14.41 The existence of the recursion can be proved in a way similar to Theo- 
rem 14.31 And the method of the new algorithm is also the same. Basically we again rewrite the left-hand 
side of the equations into powers of K^, compare their coefficients on both sides, and solve the resulting 
linear equations. Details are omitted. □ 

5. Proof of the universal denominator algorithm 

In this section we state and prove Theorem 15.11 which determines the denominator and partially the 
numerator of the rational function bj-i in Equation (jl4p . This is crucial for the performance of the algorithm 
as a whole, because it reduces the number of variables and number of equations in the final system of linear 
equations to be solved. The most straight-forward guess for the denominator 6,/_i in Equation (jl4p . i.e., 
the denominator of the right-hand side of the equation, will give us an algorithm whose performance is 
compatible to that of Sister Celine's method on a single step. Theorem IS . ll also improves Abramov's universal 
denominator [Abj . 

Let K denote a field, which for our applications it will be the field of rational functions with rational 
coefficients in a finite set of variables. Let s denote a fixed variable that does not appear in K . As usual, if 
p(s), q{s) £ K[s] are polynomials, then we write p{s) \q{s) if p{s) divides q{s). 

Consider the equation 

(15) 

where ai{s),hi{s)^c{s) G K[s\ are polynomials, and gcd(ai,6i) = 1. Define 

(t(s) — 1cm (6i(s) |0 < i < m) , 

'^ds) = Tl Tr(T(s - i),0 < I < m, 

Oi(S — l) 



and 



f(s) = gcd (to (s), ... ,T„i(s)) 

r. 

Tr{s) 



,(Pl\S + Ji) 



f{s) 



f{s) ' 



(16) Q{s) = f(s) 0,(5-1- J ), where (/..(s) 

where each Ji is the maximum of such numbers for the function 
and the outer product is over all such 0^, 



(17) R{s) = 1cm, 



0<i<rn 



h,{s - ^)Q(s) 



gcd {hi{s - i)Q{s),a^{s - i) JJ^^^^ bj(s - i)Q{s +j- i)^ 
Obviously, R{s),Q{s) G K[s] are polynomials. 

Theorem 5.1. With the above conventions, every rational solution of (jlSp has the form 

Q{s) 

where P{s) G K[s] is a polynomial. 

Proof. Suppose x{s) = with gcd{ A{s), B{s)) = 1. Then 

ai{s)a{s) A{s + i) 



So 



EUiys oys ^ys -|- tj r \ r \ 



Eai{s) A{s + i)cr(s) /X.N 
, .X TT^ = c{s)a{s), 
B{s + i) b,{s) 
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,is)a{s) n B{s + j) = E + *)^' + n + J') 

J=0 i=0 j^i 



Since Ti{s) are polynomials for all i, it follows that 



B{s + i) 

Since gcd{ A{s), B{s)) = 1, it follows that 



A{s + i)n{s + i)Y[B{s + j) 



B{s) 



n{s)Y[B{s+j ~i) 



Write B{s) = Hiio Tljlo M^+j) IlY=o9j{s), where U,V^,W are constants; and gcd(5i(s),gj(s + i)) = 1 for 
any i,j,L; and gcd{gj{s),fi{s + l)) = 1 for any i, j, and — m < I < m. We call the functions (7j(s) singletons, 
and {/i(s + j)}o<i<v, chains, in which /i(s) are the heads of chains, and /i(s + Vi) the fai/s of chains. So 
we are writing B[s) uniquely as a product of chains and singletons. 
There are two cases: 

Case I: The tail of one chains is always far apart from the head of another in B{s), i.e., gcd(/i(s + 
Vi), fjis + v)) = 1 for all 0<i,j <U and < v < m. Then 



u w 

i=0 j=0 

u w 

1=0 j=0 

w 
j=o 



Bis) 



gcd{Bis),UT=iB{s+j)) 
B{s) 



gcd{B{s),UT=iB{s-j)) 
B{s) 



gcd{B{s),l\j^.^B{s + j-t)) 
Thus the singletons have the property 



gcd(S(s),ro(s)), 
gcd(S(s),r,„(s)), 
gcd{B{s),Ti{s)),i ^ 0,m. 



^cd {Ti,0 < i < m) 



At the same time, the heads of the chains fi{s) in B{s) are factors of tq, and the tails fi{s + Vi) factors 
of Tr- Therefore each chain in B{s) factors YljLo + .?) some I. Recalling the definition of Q{s) from 
Equation (fT6|) . it follows that B{s) divides Q{s). 

Case II: The heads and tails of chains are close, i.e., gcd(/j(s + Vi), fj{s + v)) ^ 1 for some < i, j < U 
and —TO < ti < TO. In this case, 0^=0 0^=0 '^'('^ + j) ^i^l contain a chain whose head is fi{s) and tail is 
fj{s + Vj) in Q{s). This is a longer chain than what B{s) really needs, but it still guarantees that B{s) 
divides Q{s). 

So far, this proves that x{s) = 0^ where A{s) G K[s] is a polynomial. To finish the proof, it suffices to 
show that R{s) (given by Equation (|17p ) divides A{s). Since 



E ai[s) A{s + I) _ 
his) Q{s + ^)~''^''' 
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with ai,bi,c polynomials, any polynomial factor that appears only once in the m + 1 denominators on the 
left-hand side must also divide the corresponding numerator, which means 



his)Q{s + i) 



gcd [b^{s)Q{s + i), a,{s) U,^^ b,{s) Q{s + j) 



A{s + i). 



□ 



Remark 5.1. When m = 1. Theorem 15.11 becomes Gosper's algorithm. Recall that Gosper's algorithm tries 
to find rational solution x(s) such that 

a(s)c(s + £) , 

for some integer i. Based on our propositions, we get a chain Jli=o '^('^ ~^ *) ^ ^^'^ denominator and b{s — 1) 
as part of the numerator, which agrees with Gosper's result. 

Remark 5.2. Abramov's universal denominator treats the singletons in Theorem 15.11 as chains of length 
1, and then tries to find all chains. However, by picking singletons out first, we reduce the possibility of 
generating redundant chains in the denominator, because factors in the leading coefficient may mingle with 
the singletons and generate unwanted factors in chains. We illustrate the effect by example. 



Example 5.3. This is Example 1 in [Ab| . 

{n + 4)(2n + 1)(?7, + 2)x{n + 3) - (2n + 3)(n + 3)(n + l)x{n + 2) + n{n + 2){2n - 3)x{n + 1) 

-{71 - l){2n - l){n + l)x{n) = 0. 

Abramov's algorithm gives the denominator u{n) =^ ~ n for all rational function solutions x{n) S Q(n) 
of the above equation, and computes the general polynomial solution C(2n^ — 3n). However our algorithm 
finds two singletons (n-|- l)(n— 1) and no chains. So the denominator is Q{n) ~ n? — which strictly divides 
uin). 



Example 5.4. In one of the intermediate steps for Example 13.31 we get 

(fci + 1 - n + k2f{n + fc2 + 3)(n + fc2 + 2)(-n + fc2 + 2)(-?i + fc2 + 1) 



. , , -(2fci + fc2 +4A:2fci - efen - 3n + fci +3?!^ - 6fcin + 2A:?)(7i + fc2 + 2)(-7i + fc2 + 1) 

^^"''^''^'^ + (fc2 + 2)2(fci + l-n + fc2)(fci-3n + fo) x[nM,ko_ + l) 



(k2 + 3)2(fc2 + 2)2(A:i + 2-n + fc2)(fci - 3n + ^2 + 1) 

c(n, fci, fe) 



x{n, ki,k2 + 2) 



{n + k2 + l)(-n + fca) n,=o K^i - 3n + fc2 + j){ki -7i + k2 + l+ j){ki + 1 + j] 
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to solve for x{n,ki,k2) with c{n,ki,k2) a polynomial. After multiplying both sides by the denominator of 
the right-hand side, we find four singletons (n + ^2 + l){—n + k2){ki — 3n -I- k2){ki — rt + A;2 + 1); no chain 
in the denominator of x{n^ fci, fc2); and /c|(fc2 + 1)^ factors of the numerator of x(n, ki, k2). Hence 

^''^'^ = (n + + l){-n + k^)tt-in^ k2){k, n + k^ + 1) 
where P{n^ fci, ^2) is a polynomial. 

Our method keeps finding the best possible denominators in all the steps of the examples discussed in the 
paper. 

References 

[Ab] S.A. Abramov, Rational solutions of linear difference and q-difference equations with polynomial coefficients, Proc. 

ISSAC'95, ACM Press, 1995, 285-289. 
[AZ] M. Apagodu and D. Zeilberger, Multi-Variable Zeilberger and Almkvist-Zeilberger Algorithms and the Sharpening of 

Wilf-Zeilberger Theory, Adv. Appl. Math 37 (2006) 139-152. 



12 



STAVROS GAROUFALIDIS AND XINYU SUN 



[Fas] Sister Mary Celiner Fasenmyer, Some generalized hypergeometric polynomials, Ph.D. dissertation, University of Michi- 
gan, November, 1945. 

[GSl] S. GaroufaUdis and X. Sun, The non- commutative A-polynomial of twist knots, preprint 2008 l arXlv : 0802 ■ 40741 

[GS2] and , Computer data available at www.niath.gatech.edu/~stavros/publlcatlons.html. 

[Gos] R. W. Gosper, Decision procedures for indefinite hypergeometric summation, Proc. Nat. Acad. Sci, 75 (1978), 40-42. 
[PSc] P. Paule and C. Schneider, Creative telescoping for hypergeometric double sums, Preprint, 2005. 

[PSh] P. Paule and M. Schorn, A Mathematica version of Zeilberger's algorithm for proving binomial coefficient identities, 

J. Symbolic Comput. 20 (1995) 673-698. 
[PWZ] M. Petkovsek, H.S. Wilf and D. Zeilberger, A = B, A.K. Peters, Ltd., Wellesley, MA 1996. 
[Sch] C. Schneider, A new Sigma approach to multi-summation. Adv. Appl. Math 34 (2005) 740-767. 

[Weg] K. Wegschaider, Computer generated proofs of binomial multi-sum identities. Diploma thesis, RISC Linz, Johannes 
Kepler University, May 1997. 

[WZ] H. Wilf and D. Zeilberger, An algorithmic proof theory for hypergeometric (ordinary and q) multisum/integral iden- 
tities, Invcntiones Math. 108 (1992) 575-633. 
[Ze] D. Zeilberger, Sister Celine's technique and its generalizations, J. Math. Anal. Appl. 85 (1982) 114-211. 

School of Mathematics, Georgia Institute of Technology, Atlanta, GA 30332-0160, USA, |http: //www .math. gatech| 
. edu/~stavros 

E-mail address: stavrosOmath. gatech.edu 

Department of Mathematics, Tulane University, 6823 St. Charles Ave, New Orleans, LA 70118, USA, [ht"tp://| 
www . math . tulane . edu/ ~xsunl 

E-mail address: xsunlOtulane . edu 



